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Abstract 

Suppose you and your friend both do n tosses of an unfair coin with probability of 
heads equal to a. What is the behavior of the probability that you obtain at least d 
more heads than your friend if you make r additional tosses? We obtain asymptotic 
and monotonicity /convexity properties for this competing probability as a function of 
n, and demonstrate surprising phase transition phenomenons as parameters d, r and a 
vary. Our main tools are integral representations based on Fourier analysis. 
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1 Introduction 

Suppose you and your friend both do n tosses of an unfair coin with probability of heads 
equal to a. What is the behavior of the competing probability that you obtain at least d 
more heads than your friend if you make r additional tosses? 

For a fair coin with a = 1/2 and r = d = 1, this is Example 3.33 on page 118 of the 
textbook pj : "Adam tosses a fair coin n + 1 times, Andrew tosses the same coin n times. 
What is the probability that Adam gets more heads than Andrew?" 

Two solutions are offered in the textbook. One uses the symmetry and finds the answer 
1/2 easily. The other is a direct computation, following the identity 

n n+1 / , 1 \ / 

n + 1 \ (n 
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and the author noted that "a combinatorial solution to this problem is neither elegant nor 
easy to handle" . 

The same problem also appeared as Problem 21 of the self-test problems and exercises 
on page 115 of the textbook [7], and asks the following: "If A flips n + 1 and B flips n fair 
coins, show that the probability that A gets more heads than B is 1/2." This problem is at 
the end of the chapter on conditional probability and independence, and with the hint that 
one should "condition on which player has more heads after each has flipped n coins." 

What happens if the coin is unfair is not mentioned in both textbooks. In this paper, 
we consider competing probability of two independent binomials and the associated phase 
transition behaviors as parameters vary. 

To pose the problem formally, let Sk and S' n be independent binomial random variables, 
that is, Sk = Xi + ■ ■ ■ + X k and S' n = X[ + • • • + X' n , where {X i} A]}j>! are independent iden- 
tically distributed random variables that equal 1 with probability a and with probability 
1 — a. The textbooks consider the probabilities 

p n :=F{S n+1 >S' n + l} 

while we are interested in the more general 

p r /:=¥{S n+r >S' n + d}, 

where of course p^ 1 = p n ; we always assume that a 6 (0, 1) and r,d> 1. 
Clearly we can write a combinatorial expression 

n n+r 

fn = EE P ^ = 3) ■ ^(S'n = i) 
i=0 j=i+d 

= £ £ (T)(l) ait,(1 "" r_H (2) 

i=0 j=i+d ^ J / \ / 

by using independence and binomial probabilities. Unlike (pQ), the expression (j2j) for p^ d is 
of a little use in analyzing finer behaviors as parameters vary. On the other hand, by the 
Cental Limit Theorem, for any r, d and a as above, 

lim p r > d = lim > ^L±^\ = p{iV > N'\ = -, (3) 

n^oo n->oo I yn y/fl J 2 

where N and N' are independent normal random variables with mean and variance a(l— a). 
This standard technique allows one to find the limit but tells nothing about the mode and 
exact rate of convergence. Our first result provides a useful integral representation for p^f 
which implies precise asymptotic. 

Theorem 1. For any positive integers r and d, 

f/ = ~ + - T - X 2 )-V 2 P r /(x)dx (4) 

2 7T Jo 
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where 
and 

As a consequence, 



Q a {x) := 1 -4a(l -a)(l - x 2 ) (5) 



^)==y:^(^(i -«+«*)-) . ( 6 ) 



3=0 



(2j)! 



hm K' ~7i)= a r- ^ — (7) 
^ 4V7T y/a(l — a) 

The actual degree of P^ ,d (x) is 2max(r — d, d — 1) as the higher coefficients in (jHJ) 
vanish. In Section [2] we give another convenient formula for P£ d (x) in terms of Chebyshev 
polynomials of the second kind, and then easily show that 

P r /{1) = 2ar-2d + l. 

Note that for any a G (0, 1), the function Q a (x) is increasing on [0, 1] with Q a (Q) > and 
Q a (l) = 1. Therefore as n increases, the main contribution to the integral in (jlj) comes from 
a decreasing small neighborhood of 1, and a standard analysis of (|1J) implies (|7J). Another 
advantage of the integral representation (Tjf]) stems from the isolation of the variable n from 
parameters r and g?. 

Since the integral representation for p^ d in (pEJ) can be integrated out via trigonometric 
substitution a; = cos(t/2), it seems possible to check the equivalence of (j2J) and (j4]) by pure 
algebraic manipulations. However, such an approach provides no probabilistic insights into 
the integral representation and no clues on how we discovered it. We will present a proof 
via a combination of one-step-back analysis and Fourier analytic methods, along the line we 
initially derived the representation. 

Next we consider the mode of convergence in terms of monotonicity and convexity 
properties in the simplest setting r = d = 1. Let us agree that "increasing/decreasing" 
stand for "strictly increasing/ strictly decreasing" throughout this paper. 

Theorem 2. The sequence p n = p 1 / = F{S n+1 > S' n + 1} = F{S n+1 > S'J, n = 1, 2, ... , 
is monotone and convex/ concave. Precisely, p n is increasing and concave when a < 1/2, 
decreasing and convex when a > 1/2, and equal to 1/2 for all n when a = 1/2. 

It is remarkable that the sequence p]{ 1 is monotone starting from its first term. We 
first discovered this feature when studying some very specific probabilistic properties of the 
so-called double-sided exponential random variables, see the discussion at the end of the 
paper. 

Observe that as a decreases from 1 to 0, the behavior of the sequence p^ 1 changes 
instantly at the critical value 1/2, see Fig. [TJ This phenomenon becomes much clear in the 
setting with r > 2 and d — 1. In order to give the best statement of our results, from this 
point on we assume that n > rather than n > 1, which was natural for the introduction 
of the problem. We say that a sequence a n ,n > 0, is unimodal with the mode N > 1 if 

QjQ < Qi < • • • < flAT— 1 ^ CLn > > > • • • • 
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a = 0.3 



a = 0.5 



a = 0.6 



0.53 



0.5 



0.44 



Figure 1: The sequence p n , < n < 50, for different a. The dashed line corresponds to 0.5. 
Theorem 3. For r > 2, the sequence p^ 1 is increasing when < a < unimodal when 



We observe the peculiar "phase transition" in the behavior of p"^ 1 with decrease of a, 
see Fig. [2j For large a the sequence is decreasing while as a gets smaller reaching the critical 
value l/(r + 1), p 7 ^ 1 becomes unimodal and the point of maximum drifts to the right as a 
decays to l/(2r), as shown in Fig. [2cl When a reaches the critical value l/(2r), p 7 ^ 1 becomes 
increasing, which corresponds to the limit case N a ^ r — oo. It is indeed remarkable that in 
the transitional mode the sequence p 7 ^ 1 is unimodal. Relations (IE|) and fl9]) of Theorem [3] 
describe the "speed" of phase transition near the critical value a = l/(2r). 

For the general setting as parameters a, d, r vary, we have the following slightly less 
precise results which still capture the phase transition phenomenon between the monotone 
modes of convergence. 

Theorem 4. Depending on the relations between r and d, the sequence p^ d is 

1. increasing (for any 0<a<l)ifl<r<d— 1; 

2. (a) increasing (for all n) when < a < d/(r + 1) ; 

(b) increasing for large n when d/(r + l) < a < (2d — l)/(2r), 

(c) decreasing for large n when (2d — 1)/ (2r) < a < 1, 
if d < r < 2d - 2; 

3. (a) increasing when < a < 1/2, 

(b) identically equals 1/2 when a = 1/2, 
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while for the maximum 




as a 



4 



= 0.3 



a- = 0.1 



0.64 



0.56 



0.46 



0.38 



(a) 



(b) 



a = 0.214 



0.53 



0.507 - 



a = 0.184 



a = 0.174 



0.502 



0.52- 



0.503 



0.499 



(c) 

Figure 2: The sequence p^ 1 , < n < 50: monotone modes (a) and (b) and the phase 
transition (c) between them with decrease of a. In the transitional mode (c) the point of 
maximum drifts to the right. 



(c) decreasing when 1/2 < a < 1, 
ifr = 2d-l. 

4- (a) increasing for large n when < a < {2d — l)/(2r), 

(b) decreasing for large n when {2d — l)/(2r) < a < d/(r + 1), 

(c) decreasing (for all n) when d/{r + 1) < a < 1, 
ifr> 2d. 



a = 0.75 0- = 69 a = " ^ 

0.6 \ 0.45- 

0.499 



0.54h 

.49 4 | •.........-•••"* 037 •; 

Figure 3: Phase transition in the behavior of p^ 4 , < n < 50, with decrease of a. 



Part 1 is in some sense degenerate while Part 2 uncovers another type of behavior of 
p^ d , see Fig. [31 However, the reader may concentrate on the case r > 2d— 1, which is already 
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familiar from Theorems |2] and [3j as all the results Part 2 for d < r < 2d — 2 follow from Part 
4 by a certain duality relation given in Section [3j 

2 One-step-back analysis and Fourier method 

It is natural to look into the increments of p^ d in n. We start with some elementary but 
very useful consideration, one-step-back analysis. 

Lemma 1. For any r,d > 1, 

i=Q 

where the sequence, defined by 

qf :=F{S' n -S n = i}, 

is decreasing in i for % > and satisfies qn = q$ . 
Proof. We write 

n 

P r n d = F{s n+r >s' n + d} = f{s'; >s' n -s n + d}=Y, q$Hs r >i + d) 

i=—n 

with S" := S n+r — S r and similarly, 

n 

= p{s? > + -s n - x n+l + d} = J2 q%H s r >x[-x>; + i + d}. 

i=—n 

Hence 

n 

Vn d +i-Vn d = a{l -a) £ q { nOP{Sr > 1 + d + l) + ¥{S r > i + d- 1} -2F{S r > i + d]\ 

i=—n 
n 

= a(l-a) £ ^(W'SV = « + - 1} -P{5 r = i + d}), 

i=— n 

which yields the first statement of the lemma. 

The relation q^ = qn is obvious. Let us use induction to show that q$ > qn +1 ^ for 
any n and < i < n. The initial case n = 1 is trivial. We use the law of total probability 
to get from n to n + 1: 

& = «a - «) (?r x) + + (« 2 + (i - «) 2 )e 
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Then 

& - & = «a - «) (?r 1} - ^ +2) ) + (« 2 + (i - «) 2 - «a - «)) - 

which is positive by the induction hypothesis. 

□ 

In order to gain more quantitative information, we employ the powerful Fourier method. 
Recall that for any integer- valued random variable Z, it holds that 

F{Z = 0} = -!- J Ee ltz dt. (10) 

Indeed, 

Ee itz dt = e itk F{Z = k}dt = / e itk dt¥{Z = k} = 2n¥{Z = 0}. 

Lemma 2. For any r,d > 1, 

r <f-l 

:= = O^CM)^) " E P {^ = *}^(d-i-i)(a;) (11) 

is defined in terms of Chebyshev polynomials Uk(x) of the second kind. 

Recall that Chebyshev polynomials of the second kind, defined by £4(x) = *E^±pi with 
x = cos 9, are orthogonal on [—1,1] with weight y/i — x 2 . It is remarkable that P£ d {x) 
is explicitly expressed in terms of C4(x), whose properties are well known. In the next 
section we will show that F£ d (x) = P£ d (x), and thus ( 1TTT) actually serves as another useful 
representation of the polynomials defined in Theorem [TJ 

Proof. With the reminder q { n k) = ¥{S' n - S n = k}, we use (HOP to write 

= _L [* Ee u{s "- s '"- k) dt = ^ [ W ^(t)e~ lkt dt, 

where 



2tt /_„ 2vr 



yj a (t) : = (l - a + ae u ) (l - a + ae~ u ) 

= 1 - 4a(l - a) sin 2 - (12) 



is the characteristic function of S\ — S[. By the symmetry of <p a {t), 



71 

and combining ( I13p with Lemma [TJ we get 



rf° = - / cos ktdt. (13) 



- p** = a ^ a ) r ^( t )fp{5 r = z} (cos(2 -d + l)t- cos(z - d)t)dt. (14) 
71 i=0 



As for > it holds that 



., , , 2 t sin(fc + l/2)i 2 t ft 

cos(k + l)t — cos kt = —2 sin - ■ . — = —2 sin - • (7 2 fc ( cos - 

Zi siri \l I z \ _ \ z 

T4"l) transforms to 



P& - I*' = -^-^ f sin 2 1 ■ i^(cos l)dt, 

J 

and we conclude the proof with the change x = cos(t/2). □ 



3 Proofs of the theorems. 

Proof of Theorem [0 First check that 

P?{x) = P£*(x). (15) 

Using the representation 

u 2k ( X ) = YX-iy ( 2k ~ j ) w 2k - 23 = E(- 1 ) fe " ( k t j ) (2x)2j 

j=o \ 3 J j=0 V 47 / 

from Section 6.10.7 of [8], we obtain 

p r a d w = EE p {^ = 4(-r- d ( l+ [^)(2x) 2j 

i=d j=0 \ 3 / 

- E = oc- 1 )^- 1 ( d ~ * 2 + j " l ) 

i=o j=o \ 3 / 

Let us agree that := ^ rim=n-fc+i m anc ^ (o) := 1 to wri te binomial coefficients with 
negative n. This allows us to note that ( -t ^' -1 ) = ( I+ 2~ ) f° r an y d > 0. The next 
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important observation is that in both double sums the summation in j could be taken from 

to oo. Indeed, any j > i — d + 1 gives no contribution to the first sum as the product 
rim=^-i-d+i m ' wri ich corresponds to the binomial coefficient, includes a zero factor because 

1 + j — d > and i — j — d + 1 < 0. The same applies to the second double sum for j > d — i 
as the product Yltn=d-7-j 171 * s zero smce d — i + j — 1>0 while d — i — j < 0. Thus 



r oo / * _|_ ' _ rl\ 

i=0 j=0 \ J / 

^ dt 2 i V 



2.7 



(2J)J 



Then (|15p follows since the generating function of 5* r is Et 71 = (1 — a + crf) r 
It now remains to use Lemma [2] and (fTTj) to write the telescoping sum 



i=0 " ^° i=0 

r.d 



As fc — > oo, T^n+k ~~ ^ 1/2 by ([3]), and simplifying the sum of the geometric series with (j5]), we 
get (jlj) by the dominated convergence theorem. 

For the precise asymptotic relation (j7|), we return to the convenient trigonometric sub- 
stitution x = cos(t/2) in (jlj) and get 



"5 = sf rtW^O"!)* 



with yj a (t) = Qa(cos(t/2)) = 1 - 4a(l - a)sin 2 (t/2). Choose a 5 > such that t/4 < 
sin(£/2) < t/2 on [0,5] and observe that <p a (t) is decreasing on [0, tc] and <p a (0) = 1. We 
have 



^(t)^(coB|)* = / 5 ^(t)P Q M (cos^)rft + 0«(5)) 

«/ 

= ^jT^M"^)** ^ (16) 

lim 99" ( -4= ) = lim ( 1 - 4a(l - a) sin 2 — ^= ] 



where 



Now 



e -a(l-a)^_ ^ 



lim ri [0 ^( s )^(^)p Q ^(cos-^)rf S = ^ e~< 1 -^ 2 r/(l)ds 



7T 



2 V a(l - a) 



P r, d{1) 



by the dominated convergence theorem, with e a<yl a ) s2 / 4 -sup 0<:r<1 \P^ d (x)\ as an integrable 
majorant. Combining the arguments above together, we finish the proof of Theorem 1 once 
we find P r /{1). 

It is well known that £4(1) = k + 1, hence (fTTT) and (f!5l) imply 



P£<*(1) = ^P{ 1 S r = i}(2z-2d+l)-^P{5 r = z}(2c/-2z-l) 

i=d i=0 

= E(2S r - 2d + 1) = 2ar - 2d + 1. (19) 

□ 

Proof of Theorem El As P^' 1 (x) = 2a — 1, Lemma [2] yields 



. - Pn = 4a(1 - a)(1 - 2a> /" <K(x)^ato. 
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It is readily seen that p n is monotone because the integrand is nonnegative; moveover, the 
integrand is monotone in n for each x, so p n +i ~Pn is monotone implying convexity of p n . □ 

Note that monotonicity of p n could be obtained directly from Lemma [H whose proof 
requires only an elementary one-step-back analysis. The same is true for convexity of p n but 
some additional study of properties of q$ should be done. 

It is also worth mentioning that the asymptotic of p n — 1/2 could be found via purely 
probabilistic argument with no use of Fourier method. Indeed, the trivial identities 1 = 
2¥{S n > S n } + ¥{S n = S n } and ¥{S n+1 > S n } = a¥{S n = S n } + F{S n > S n } imply 

Pn = \ + ^Y^V{S n = S n }. (20) 

The asymptotic of the probability in the right-hand side is given by the classical local limit 
theorem. 

We now turn to the proof of Theorem H] as it partly covers Theorem 

Proof of Theorem 31 In the case 1 < r < d — 1, the result immediately follows from 
Lemma [TJ In the case of r > 2d — 1, we start with the proof of Cases (3c) and (4c). By 
Lemma [H we write 

,d-l 

Pnti-P r / = a(l-a)^(F{S r = z}-F{S r = 2d- i -l}){qt d+1) -qt d) ) 

+ irr{s r = i }( q t d+1) - q t d) ))- 
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The second sum is non-positive. Let us show that 

¥{S r =2d-i- 1} 



¥{S r = i} 



> 1, < i < d- 1 



when a > d/(r + 1) to prove that the first sum is negative. This statement generalizes the 
well known result that the (last) maximum of binomial coefficients P{SV = z} occurs at 
i = [(r + l)a]. 

Consider the ratio 

a i+1 _ (1 - a) 2 (z + l)(2d - % - 1) 
aj a 2 (r - i)(r - 2d + i + 2) 

and rewrite it in the form 

a\2 d 2 -{d-i-l) 2 /l-a\2/ ( r _d + i)2_ rf 2 



1 



a ) (r -d+ l) 2 - (d-i - l) 2 V a / V (r - d + l) 2 - [d - i - 1] 

to observe that this quantity increases in i for < i < d — 1 when r > 2d and is constant 
when r = 2d — 1. Now as a > d/(r + l), in both cases we have 

a i+ i < a d _ / l-« \ 2 d 2 ^ i 



a, a^-i V a / (r — d+1) 2 
and thus Oj is decreasing. Then 

F{S r = d\ 

a i > a d-i = — r - ^ > 1 

P{S r = d- 1} 

because a; > d/(r + 1). 

Next we present the proof of Cases (3a), (4a), and (4b). First note that a < d/{r + 
1) < 1/2. Then by Lemma [21 the sign of pjf+i — Pn d for large n is opposite to the sign of 
P^(l) = P£ d (l) = 2ar -2d + l, see and P]). 

The case a = (2d — l)/(2r) requires more attention. Here P^ d (l) = 0, and the sign of 
Pn+i—Pn d f° r l ar g e n coincides with that of -j^P£ d (l). We use the formulae T^(x) = nU n -.\(x) 
and T"(l) — (n — l)n 2 (n + l)/3, where T n (x) = cos(n arccosx) are Chebyshev polynomials 
of the first kind, to find that U' n {l) = n(n + l)(n + 2)/3. Hence, arguing as in the proof of 
(|T9|) . we obtain 

^-P r ' d (l) = -E(2S r - 2d)(2S r -2d+ l)(2S r - 2d + 2). 
ox 3 

As the expectation in the right-hand side equals 



Q3 



t=l 



t=l 
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we substitute a = (2d — l)/(2r) and after some simplifications get 

d_ d _ 2(2d -!)(! + + 3r + 2r 2 - 2d(2 + 3r)) 

^■'(ad-ij/carjV 1 ; - 3r 2 ■ ^ 

In order to check that this expression is positive, note that the partial derivative in r of the 
last factor in the numerator is equal to 4r + 3 — Qd > as r > 2d. Hence the minimal value 
of this factor is attained at r = 2d and equals 2c? + 1 > 0. 

In the case of d < r < 2d — 1, we observe the duality relation 

p r /(a) = l-p r /- d+1 (l-a), (22) 

which follows from 

p r /(a) = F{S n+r (a) > S' n (a) + d] = F{n + r - S n+r (a) <n- S' n (a) + r - d} 

by comparing tails. Here we temporarily changed the notation to stress that p^ d is a function 
of a. 

Now Case (3b) follows immediately while (122]) implies that Case (2) is equivalent to Case 
(4), which was proved above. Indeed, if d < r < 2d — 2, then r > 2d', where d' :— r — d+ 1; 
and if r > 2d, then d' < r < 2d' — 2. Thus the function (r, d) t- > (r, r — d + 1) is a bijection 
between the sets {(r,d) : d < r < 2d-2} and {(r,d) : r > 2d} while l-d/(r + l) = d'/(r + l) 
and 1 - (2d- l)/(2r) = (2d'- l)/(2r). 

□ 

Proof of Theorem [3]. It is well known that the maximum of £4(x) on [0,1] occurs at 
ac = 1. Hence (fTT]h ([15]) and (fig]) imply that P^(x) < P r /{1) = 2ar - 1 for x G [0, 1], and 
by Lemma [21 the sequence p 7 ^ 1 is increasing when < a < l/(2r). The decrease of p^ 1 for 
l/(r + 1) < a < 1 is already covered by Case (4c) of Theorem |H 

For a proof of unimodality in the transitional zone it suffices to check that P^ ,l (x) has 
only one root on [0,1] when l/(2r) < a < l/(r + 1). Indeed, let P^ 1 (xq) = for some 
x G (0, 1), and let — p 1 ^ 1 < for some k > 0. We claim that for any n > k it holds that 

Q n Q (x)VT^P r a d (x) > Q n a - k (x )Q k a (x)VT^P^ d (x), x G [0, 1] \ {x }. (23) 

First observe from © that P^{0) = -(1 - 2a) r < while P r /(l) = 2ar - 1 > 0, so 
P£ l (x) < on (0,xo) and P£ l (x) > on (xq,1). Then we get (1231) using that Q a (x) is 
positive and increasing in x on [0, 1]. Now integrate ( 1231 over [0, 1] and apply Lemma [2] to 
get Pn+i —Pn 1 < Qa~ k { x o){Pk+i ~ Pk ) — f° r a ^ n > k. Thus unimodality follows from the 
uniqueness of root. Note that by Lemma (H 

- p f = a[l -a)(F{S r = 0} -F{S r = l}) > 

assuring the the mode satisfies N a ^ r > 1. 
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Clearly, we prove that P^ ,l (x) has only one root on [0, 1] when l/(2r) < a < l/(r + 1) 
if we check that P^ 1 {x) is increasing on [0, 1] when a < l/(r + 1). We claim the stronger 
statement: with exception of the constant term, all coefficients of the polynomial P^ 1 {x) are 
positive when a < l/(r + 1). 



By Theorem (TJ we should show that *(1 — a + at) r ) 

1 < j < r — 1 . The Leibnitz formula gives 



i=-l 



is positive for any 



dt 2 i 



2jAr 

(t j -\l-a + at) r ) = (- 1 

k=j+l 



\k-j-l 



a k, 



where a k := ( 2 ^) ( r _^i| fc ^ ! _ 1 ) ! a fc (l — 2a) r fc , and it suffices to prove that are decreasing. 
We have 

Qfc+i a (r-k)(2j-k) 
a k ~ I -2a (k-j)(k + l) ' 
which is obviously decreasing in k, so 



a k+l ^ O-j+2 



a (r-j — l)(i — 1) 



2a 



J + 2 



< 



a 



2a 



r - 2 
[r - 2) < < 1 



when a < l/(r + 1). This completes the proof of unimodality of p^ 1 . 

To prove ([H]), we use the argument similar to the proof of ([7]) in Theorem [TJ By Lemma [2] 
and the change x = cos(t/2), we have 



<(t)sin 2 ^P^(cos^)dt, 



2a(l - a) 

and representing P^ 1 {x) via its Taylor polynomial at x — 1 



2a(l - a) 



ipl(t)sm 2 -P r /(l)dt 



+ 



^(t)sin : 



^ 1 < i )( i -»4 



i d 2 



2dx 2 



P^(9(t)){l-cos 



t\ 2 
2 



dt 



--: -h(a,n) + I 2 (a,n) 



where cos(t/2) < 6{t) < 1. 

Now assume n = n(a) — > oo as a I l/(2r). We claim that 



Ji(a, n(a)) 



2ar - 1 



2(2r - l) 3 /2 n (a) 3 / 2 



y^F(r - l)r 3 1 1 
J 2 (a,n(a))~ 4(2r _ i)3/2 --^ as a|-. 

(24) 
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Indeed, arguing as in ( !T6|) and ({TBI) and using (|T7|) . which of course holds when a —¥ l/(2r), 
we get 

»(«» ~ y o ~ i ^ 

by the dominated convergence theorem. It only remains to compute the integral and use the 
equation P ( ^' 1 (l) = 2ar — 1 to get the first part of ff24l) . Similarly, 

T I I \\ fe-^Q ^> / n( S \ dx^l/(2r)y l ) / -2l=I s 2 4 

^(«.«(«)) - -^r i d S ~ 32n5/2 y o e ^ S rf S 

as the term coming from the second derivative of P'^ 1 (x) gives a lower order contribution, 
and we get the second relation in (1241) by computing the integral and using (j2~T|) . 

Clearly, A^, )r -> oo as « j l/(2r), and ( 124"1) implies Ik(a,N air — 1) ~ N ajr ), k = 
1,2. By the definition of iV air , we have l2(a,N air — 1) > Ii(a,N ajr — 1) and Ii(a,N ajr ) > 
J 2 («, N a>r ), hence 

7i(a, A^ a>r ) ~ J 2 (a, A^,r) as a 

2r 

We now apply (|24|) again to get (J8j. 

For ([nD, we argue in the same manner as above. Omitting the details, we get 

and the required relation follows. □ 



4 Additional remarks and related questions 

Here we mention some issues that deserve more attention than we can provide. Then we 
explain our initial interest to the problem. 

Simulations show that as in the case p 7 ^ 1 with r > 2, the sequence p^ d with r > 2c? is 
increasing when < a < (2d— l)/(2r) and unimodal in the transitional phase (2d— l)/(2r) < 
a < d/(r + 1). The proofs presented here do not admit a reasonable generalization even for 
d = 2. For example, simulations show that the coefficient of P 7,2 (x) at x is negative and all 
the other coefficients are positive. This could be neither easily proved nor used for a proof 
similar to the one given here for d = 1 as P 7,2 (x) has two roots on [0, 1] for some a. Once 
unimodality in the transitional phase is established for general d, the asymptotic of the mode 
N£ d could be found by exactly the same argument as used for (JHJ. 

Another open question is to verify that the sequence p^f is convex/concave in the mono- 
tone modes, as observed by simulations. By the same argument as in the proof of Theorem [TJ 
we easily get convexity (concavity, to be precise) of p]] 1 when < a < l/(2r). Unfortunately, 
the same method does not seem to work for l/(r + 1) < a < 1. 
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Next we should mention that, after a preliminary draft of this paper was finished, we 
found that a related problem was studied in pQ . Their analysis uses the multivariate form of 
Zeilbergers algorithm to find representations and estimates based on a linear combination of 
two Legendre polynomials. Also, after submitting this paper, Tamas Lengyel and the referee 
pointed out the paper [6] on approximating point spread distributions. Our Fourier analytic 
method can also be used to simplify some of the arguments in both [1] and [6], and provide 
integral representations for general r, d, p and q, see [2] for details. 

Finally, let us explain the origins of our interest to the problem considered in this paper. 
We say that a random variable Z has a double-sided exponential distribution if it has a 
density p(x) of the form p(x) = aae~ ax for x > and p(x) = (1 — a)be bx for x < 0, where 
a, b > and a G (0, 1). Such a specific random variable has very nice properties from the 
point of view of fluctuation theory, see [U Ch. VI. 8]. KZ = means b = ^—^-a so any centered 
double-sided exponential variable could be written in the form Z = \{XY — t^-(1 — X)Y^j, 
where X is a Bernoulli random variable that equals 1 and with probability a and 1 — a, 
respectively, Y is a standard exponential random variable with density e~ x for x > 0, and 
X and Y are independent. 

Consider independent identically distributed double-sided exponential random variables 
Z n . Initially, our goal was to show that the sequence P{Zi + ■ • • + Z n > 0} — 1/2 does not 
change its sign as n increases. This property was verified by simulations, and a rigorous proof 
was required. One of the possible approaches is to show that the sequence P{Zi + • • • + Z n > 
Oj is monotone since by the central limit theorem, its limit equals 1/2. 

Introducing independent sequences of standard exponential random variables Y n and 
Bernoulli random variables X n , we write 

a) n - k ¥{Y 1 + ■ ■ ■ + Y k > ^— (Y k+l + ...*;)} 
a )«-* P | yi + ... + Yk >a(Y 1 + ...Y n )} 

a) n - k P{U k>nr . 1 >a} + a n , 

where we used the well known (Karlin [5]) representation of uniform order statistics in terms 
of exponential random variables. 

Recall the appropriate definition. Let U\, . . . , U m be independent random variables that 
are uniformly distributed on [0,1], and arrange them in the ascending order. The k-th 
element of this new sequence is called the k-th order statistics of Ui, . . . , U m and denoted by 
V k)Tn . By definition, Ui jm < • • • < U mtTn but all the inequalities are strict with probability 
one. 
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Note that F{U k}Tn > a} =F{S m < k}, hence 

n-l , s 

¥{Z l + • • • + Z n > 0} = \l) ak{l ~ «) n ~ fc P{^-i <k}+ a n = F{S' n > S n ^} = p n _ x . 

We see that monotonicity of P{Zi + • • • + Z n > 0} is exactly the monotonicity of p n given 
by Theorem [21 

Acknowledgement: The authors would like to thank Tamas Lengyel and the referee 
for pointing out the paper [S] and their insightful remarks and suggestions. 
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